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Abstract: We use the AdS/CFT correspondence to study the thermodynamics of massive 
M =2 supersymmetric hypermultiplet flavor fields coupled to = 4 supersymmetric SU (Ac) 
Yang-Mills theory, formulated on curved four-manifolds, in the limits of large Ac and large 
't Hooft coupling. The gravitational duals are probe D-branes in global thermal AdS. These 
D-branes may undergo a topology-changing transition in the bulk. The D-brane embed- 
dings near the point of the topology change exhibit a scaling symmetry. The associated 
scaling exponents can be either real- or complex- valued. Which regime applies depends on 
the dimensionality of a collapsing submanifold in the critical embedding. When the scaling 
exponents are complex- valued, a first-order transition associated with the flavor fields ap- 
pears in the dual field theory. Real scaling exponents are expected to be associated with 
a continuous transition in the dual field theory. For one example with real exponents, the 
D7-brane, we study the transition in detail. We find two field theory observables that diverge 
at the critical point, and we compute the associated critical exponents. We also present an- 
alytic and numerical evidence that the transition expresses itself in the meson spectrum as a 
non-analyticity at the critical point. We argue that the transition we study is a true phase 
transition only when the 't Hooft coupling is strictly infinite. 
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1. Introduction 

The Anti-de Sitter /Conformal Field Theory {AdS /CFT) correspondence is the conjectured 
equivalence between AA = 4 supersymmetric SU{Nc) Yang-Mills (SYM) theory and type IIB 
string theory formulated on the background spacetime AdS^ x [1~3]. Here, AdS^ is a 
"Poincare patch" of five-dimensional anti-de Sitter space, and is a five-sphere. The AdS^ x 
spacetime arises as the near-horizon geometry of Ac ^ oo coincident D3-branes. Taking 
large Ac, small string coupling gg, and fixed but large QsNc, the low-energy dynamics of 
the string theory is well-approximated by type IIB supergravity. In the SYM theory these 
correspond to the limits of large Ac and small Oym^ large 't Hooft coupling A = gyj^^Nc- 

Thermal equilibrium of the M = 4 SYM theory corresponds to non-extremal D3-branes, 
whose near-horizon geometry is five-dimensional AdS-Schwarzschild times [4, 5]. The 
temperature of the = 4 SYM theory coincides with the Hawking temperature of the AdS- 
Schwarzschild black hole. 

The M = 4 SYM theory only contains fields in the adjoint representation of the gauge 
group. To make the theory more closely resemble the theory of strong interactions. Quan- 
tum Chromodynamics (QCD), we may introduce massive "flavor" fields in the fundamental 
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representation of the gauge group. We will introduce a finite number Nf of them, so that 
Nf ^ Nc in the large limit, and work to leading non-trivial order in Nf/Nc- These flavor 
fields appear in the supergravity description as Nf probe D-branes [6], or "flavor branes." 
In particular, a probe D7-brane describes a massive M = 2 supersymmctric hypcrmultiplet 
propagating in 3+1 dimensions and interacting with the AA = 4 SYM fields. In the Nf <C Nc 
limit, the contribution that these flavor branes make to the stress-energy tensor is dwarfed 
by the contribution of the Nc D3-branes. To leading order in Nf/Nc, one may neglect the 
back-reaction of the flavor branes on the geometry, hence the label "probe." The flavor 
brane action is then the Dirac-Born-Infeld (DBI) action characterizing the embedding of the 
D-brane in the background spacetime. 

In the supergravity description, the embedding of the D7-brane undergoes a topology- 
changing transition as the temperature increases. This corresponds in the dual field theory 
to a first-order phase transition associated with the flavor fields [7 14]. In Refs. [11, 14], 
analysis of a general class of probe D-brane systems (based on supersymmetry-preserving 
Dp-Dq intersections) revealed that in these systems the transitions associated with probe 
D-branes are generically first order. 

In Ref. [13], two of the present authors analyzed the thermodynamics of flavor fields in 
the SYM theory formulated on a 3-sphere. The supergravity dual in this case does not arise 
as the near-horizon limit of any known D-brane construction in string theory: A/" = 4 SYM 
formulated on a 3-sphere, in the 't Hooft limit and at large 't Hooft coupling, is holographically 
dual to supergravity formulated on global AdS, as opposed to Poincare patch AdS. Prom 
the SYM theory perspective, the 3-sphere introduces a scale into the CFT allowing for a 
thermal phase transition. (Without such a scale to set a transition temperature, there can 
be no phase transition at any non-zero temperature.) At high temperature, relative to the 
inverse 3-sphere radius, the M = 4 SYM theory is in a "deconfined" phase, in which the free 
energy scales as N^. At low temperature, however, the free energy is order one, i.e., 0{N^). 
A first-order phase transition separates the two phases. In the supergravity description, the 
transition appears as a Hawking-Page transition [4, 15] in which the high-temperature phase 
is global AdS'-Schwarzschild and the low-tcmpcrature phase is global AdS with a periodic 
time direction, which we will call global thermal AdS. 

We may introduce a probe D7-brane into this background, which now corresponds to 
adding a flavor hypermultiplet to A/" = 4 SYM defined on S^. In the high-temperature phase, 
a probe D7-branc undergoes a topology-changing transition that is essentially the same as 
that in Poincare patch AfiS'-Schwarzschild [7-14]. This, once again, represents a first-order 
transition in the flavor physics of the dual field theory [13]. In the low-temperature phase the 
probe D7-brane undergoes a different topology-changing transition. The surprising result of 
Ref. [13] was that this transition does not appear to correspond to a first-order transition in 
the SYM theory. 

Our goal in this paper is to understand the nature of this transition. To do so, we 
will provide a more general context in which to understand the D7-brane's transition. We 
will examine other probe D-branes that are supersymmetric at zero temperature (and in 
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the Poincare patch), in particular a probe D5-brane [16-18], which describes flavor fields 
propagating in 2+1 dimensions. We will also study the AA = 4 SYM theory formulated on 
certain other curved manifolds. 

As in Refs. [11, 14, 19,20], an essential tool in our analysis will be scaling symmetries of 
probe D-brane embeddings. The crucial ingredient for the arguments of Refs. [11, 14, 19,20] 
was that certain scaling exponents associated with probe D-brane embedding solutions were 
complex-valued, which implied the existence of infinitely many solutions of the embedding 
equations (and boundary conditions). Such "multi- valued" supergravity solutions imply the 
presence of discontinuities in SYM theory observables. As external parameters, such as the 
hypermultiplet mass, are varied, the physically relevant solution can jump from one branch 
to another, leading to discontinuities in generic observables. (The details of the argument 
will be reviewed below.) 

We will find that the D7-brane in global thermal AdS has real-valued exponents, and 
that this will lead to very different physical consequences. Our general analysis will reveal 
the condition that determines whether the scaling exponents are complex or real. In the D- 
brane's topology change, we can identify a "critical solution" intermediate between the two 
topologies. The dimension of the submanifold that collapses to zero volume in this critical 
solution determines whether the scaling exponents will be real or complex. 

For the D7-brane in global thermal AdS we will further characterize the transition by 
computing two observables that are especially convenient to evaluate using the supergravity 
description: the expectation value of the hypermultiplet mass operator (i.e., the supersym- 
metric completion of ipip), and the expectation value of a particular supersymmetric Polyakov 
loop correlator. We will find that both have derivatives that diverge near the critical point, 
and we will compute the associated critical exponents. We will also compute the meson spec- 
trum of the SYM theory from fluctuations of the D7-brane's worldvolume fields. We will 
find that the transition appears in the meson spectrum as a non-analyticity, a cusp, in the 
spectrum of one scalar meson. Lastly, we will argue that the transition is only a genuine 
phase transition in the flavor physics at infinite 't Hooft coupling, that is, the divergences we 
find in derivatives of the free energy only exist when A is strictly infinite. 

This paper is organized as follows. In section ^ we present in greater detail the systems 
that we will study. In section ^ we perform the scaling analysis, compute the scaling expo- 
nents, and present numerical evidence that complex exponents signal a first-order transition 
in the SYM theory. In section ^, we compute the critical exponents and meson spectrum 
for the SYM theory dual to the probe D7-brane in global thermal AdS, and argue that the 
transition will be smoothed out at finite A. We conclude with a brief discussion in section ^. 

2. The Systems in Question 

We begin with a more precise description of the supergravity systems we will study. After 
reviewing global AdS and introducing our conventions, we describe the relevant probe D- 
brane embeddings in greater detail. 
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2.1 The Background Geometries 

We are interested in AA = 4 SYM theory formulated on curved spatial sections. We work in 
the limits of large and large 't Hooft coupling, A = gyj^iNc S> 1 so that, via the AdS/CFT 
correspondence, we may compute observables in the SYM theory using supergravity. We are 
thus interested in global AdS^^ whose boundary is a curved four-manifold. Global AdS^ does 
not arise from any known string theory construction (such as a near-horizon geometry of very 
many D-branes), so we must assume that the conjectured equivalence of theories, originally 
motivated by the black D3-brane solution of type IIB string theory [1], can be extended to 
include supergravity on global AdS^. 

As we will be studying the thermodynamics of the SYM theory, we introduce a temper- 
ature T in the standard fashion by working in Euclidean signature and compactifying the 
time direction into a circle of radius Ri = l/(27rT). We will mainly be interested in the SYM 
theory formulated on a spatial 3-sphere (times time) . We denote the radius of the 3-sphere 
as i?3. We will thus be studying the SYM theory formulated on boundary spacetimes such 
as X S^. 

We will use a global thermal AdS^ metric 

ds'^ = dp^ + cosh^p dr^ + sinh^p dVtl , (2.1) 

where we have set the curvature radius of AdS^ equal to one. We will work in these units 
throughout. In these units, we convert between string theory and SYM quantities using the 
relation a'~^ = ^TrgsNc = gyj^jNc = A, where a' is the square of the string length: a' = i"^. 
In Eq. (|2.1| ), p is the radial coordinate, r is the compact Euclidean time coordinate of period 
1/T, and dO,^ is the metric of a unit 3-sphere. The center of the AdS space is at p = 0, and 
the boundary is at /> = oo. Notice in particular that the 3-sphere collapses to zero volume at 
the center of AdS. 

We will also frequently find another coordinate system convenient, a so-called Fefferman- 
Graham coordinate system [21] in which the radial coordinate is z = e~^, so that the center 
of AdS^ is at z = 1 and the boundary is at z = 0. The global ^dSs-Schwarzschild metric, 
written in this Fefferman-Graham coordinate system, is 

ds^ = ^ [dz'^ + Qij dx'dx^] , (2.2a) 

with 

Qij dx'dx^ = \i^- ^^I^H? ^{^y^ dT^ + \ Hz) dnl , (2.2b) 

J^{z) = l-2z^ + z'^/zjj , (2.2c) 

and Zfj = [1 -|- 4(7rT)'^]~^/^ is the horizon radius. Setting = 1 produces the global thermal 
AdS^ metric, Eq. ( ^.1]) , in Fefferman-Graham coordinates. Notice that in our units the radius 
of the 3-sphere at the boundary is = 1/2. 
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When discussing global thermal AdS^, we will use the p coordinate unless stated other- 
wise. When discussing ^diSs-Schwarzschild, we will use only the z coordinate. 

A Hawking-Page transition [15] connects global thermal AdS^ and global ^dS'5-Schwarz- 
schild. This is a first order transition associated with black hole condensation. The two 
metrics written above arc both solutions to Einstein's equation with constant negative cur- 
vature and an asymptotic boundary that is x 5^. The difference of the Einstein-Hilbcrt 
action for these two solutions may be interpreted as a difference in free energies (divided 
by the temperature), and thus determines which geometry is thcrmodynamically preferred. 
Thermal AdSs is preferred (has lower free energy) at low temperatures, but a first-order 
transition occurs at a non-zero critical temperature Thp = ^ (in units of the AdS radius), 
above which AdS's-Schwarzschild is preferred [4,15]. This is a topology-changing transition: 
the topology of thermal AdS^ is x while that of AdSs-Schwarzschild is x [4, 15]. 

This transition is interpreted in the boundary SYM theory as a deconfinemcnt transition 
[4]. As the SYM theory is in finite volume, wc must be careful to define what wc mean by 
the word "confinement." A sufficient definition for our purposes may be given in terms of the 
large- ATp behavior of the free energy: if the order contribution to the free energy is zero, we 
declare the theory to be in a confined phase, whereas if the order contribution is nonzero 
we declare the theory to be in a deconfined phase. In this sense, the AT = 4 theory on x is 
in a confined phase at low temperature and in a deconfined phase at high temperature, with a 
first-order transition separating the two phases. In the absence of fundamental representation 
matter fields, this criterion is equivalent to defining a deconfined phase as one in which the 
ZjVc center symmetry is spontaneously broken in the AT^, — 00 limit, which is possible even in 
finite spatial volume. 

Notice also that the only scales in the SYM theory are the temperature T and the radius 
of the 3-sphere R3, and hence the only physically meaningful quantity is the dimensionless 
product TR3. We may interpret the limit TR3 00 either as a high-temperature limit in 
fixed volume or as a large-volume limit at fixed temperature. We will be working in a fixed 
volume, so we will interpret this as a high-temperature limit. In this limit our analysis should 
agree with any known physics of the finite-temperature SYM theory in flat space. Prom the 
supergravity perspective, our analysis should agree with any known physics in Poincare-patch 
AdS's-Schwarzschild. 

In the course of our analysis we will study the SYM theory formulated on other curved 
manifolds besides x S*^. To do so, we will exploit the fact that global AdS admits various 
"slicings." To understand these, recall that at the boundary of AdS space the metric has a 
second-order pole. To extract a boundary metric from a bulk AdS metric, we must choose 
a defining function which is arbitrary save for one feature: it has a second-order zero at the 
boundary. Multiplying the bulk metric by this defining function will produce a finite boundary 
metric. Different slicings of AdS are simply coordinate reparameterizations that naturally 
suggest different defining functions. Implicitly, we will be using these "natural" defining 
functions, which give rise to boundaries with different geometries. We will be interested in 
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slicings producing the boundary spacetimes^ AdS^, AdS^ x 5^, AdS2 x 5^, 5^ x 5^, and S*^. 
We may write the AdS^ metric for these various shcings as 

ds^ = dp^ + cosh^ p ds\j^g^_^ + sinh^ p dQ,f (2-3) 

with dQf the metric of the unit /-sphere 5' and / = 0, . . . 4. When / < 2, the shcing includes 
^^AdS4 "^l^ich denotes the AdS^^i metric. We adopt a form for the AdS^^i metric that is 
identical to Eq. (2.1), except for the replacement of the 5"^ metric, dfl^, by an S*^"' metric. 



dQ2-i- When / = 3, we adopt a convention that AdSi = S^. Notice that the S*' collapses to 
zero volume at the center of AdS. 

Not all of these boundary spacetimes have finite volume, but all introduce a spatial 
curvature scale into the dual SYM theory, allowing the otherwise scale-free = 4 SYM theory 
to undergo a thermal transition. We denote the Ricci scalar of the boundary spacetime as TZ. 
For example, in x slicing, TZ = (j/R^- 

The supergravity theory is formulated on the ten-dimensional spacetime AdS^ x (or 
AdSs-Schwarzschild times S^). Convenient forms of the metric that we will use are 

dnl = de"^ + sin^ e dnl_j + cos^ e dil] , (2.4) 

where j is any integer from 1 to 4, and 9 runs from to 7r/2. 

2.2 Probe Dp-branes 

We now introduce Nf flavor fields in the SYM theory, in the Nf <C Nc limit, which, in the 
supergravity description, corresponds to introducing Dp-branes probing the above geometries. 
Let us momentarily consider a background that is Poincare-patch AdS^ times S^, correspond- 
ing to the = 4 SYM theory in flat space and at zero temperature. This geometry arises as 
the near-horizon limit of a number Nc ^ oo D3-branes [1]. 

Known supersymmetric embeddings of probe Dp-branes in this background include D7- 
branes extended along AdS^ x S"^ [6], D5-branes extended along AdS/i x [16-18], and D3- 
branes extended along AdS^ x [22] . Introducing D7-branes corresponds in the SYM theory 
to introducing = 2 supersymmetric hypermultiplets that transform in the fundamental rep- 
resentation of the gauge group and that propagate in 3+1 dimensions. Introducing D5-branes 
corresponds to introducing flavor fields that again transform in the fundamental representa- 
tion of the gauge group, but that are confined to propagate along a (2+l)-dimensionsal 
defect. Similarly, introducing D3-branes corresponds to introducing flavor fields confined to 
propagate along a (l+l)-dimensional defect. More generally, any probe Dp-brane that is not 
extended along all of the AdS^ directions will be dual to flavor fields confined to propagate 
along some defect. In each case, the flavor fields may be given a supersymmetry-preserving 
mass m. In the supergravity description, this mass m is encoded in the geometry of the 
Dp-brane in a way that we will make explicit below. 



^The choice of defining function can also determine the topology of the boundary. For example, a defining 
function with a zero at some point on the boundary may encode a collapsing cycle in the boundary spacetime. 
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In what follows, most of our attention will be focused on Dp-branes that are supersym- 
metric in Poincare-patch AdS^ times S^. In particular, we will present numerical results for 
D7-branes and D5-branes. We would like our analysis to be as general as possible, however, 
so we will not restrict the p value of the probe Dp-branes we study. 

To introduce the Dp-brane's induced metric and action, we begin in the low-temperature 
phase, where the background geometry is global thermal AdS^ x 5^ (Eqs. and ( |2.4D ). 

The Dp-branes that we will study will be extended along AdSi x , for some i and j that 
sum to p + 1. In other words, the Dp-brane will be extended along an AdSi subspace inside 
global AdS^. This AdSi subspace includes the radial and time directions as well as an S**"^ 
inside the of AdS^. The Dp-brane will also be extended along an inside the factor. 
In the low-temperature phase, we will use the other slicings of global thermal AdS^, Eq. (|2.3D , 
only when studying D7-branes, which are extended in all of the AdS^ directions (i = 5 and 
j = 3). 

The Dp-brane may move in directions orthogonal to its worldvolume. We will allow the 



Dp-brane to move in the 9 direction defined in Eq. (2.4) and, to preserve translation invariance 
in the dual field theory, we will allow 9 to vary only as a function of the radial coordinate p. 
In other words, we allow the position of the C to vary as the Dp-brane extends in the 
radial direction. From the Dp-brane worldvolume perspective, 9{p) is a scalar field. As we 
will review in detail below, this worldvolume scalar is dual to a SYM theory operator that is 
the (supersymmetric completion of) the mass operator of the flavor fields. This is how the 
mass m is encoded in the geometry of the Dp-brane. The induced Dp-brane worldvolume 
metric is 

dsj)p = [1 + 9'{pf] dp^ + coshV dr^ + sinhV dn'f_2 + cos'^9{p) dQ] (2.5) 

where the prime denotes differentiation with respect to p. In the high-temperature phase, the 
embeddings have the same form: AdS'j-Schwarzschild times C S^, and now the worldvol- 
ume scalar 9{z) depends only on z. 

To determine 9{p), we need an equation of motion and boundary conditions. The equation 
of motion comes from the probe Dp-brane action, the Dirac-Born-Infeld (DBI) action,^ 

Sdp = NfTopJ d(P+i)C \/det(5afe) , (2.6) 

where T^p = {27r)^Pg^^a'^^P'^^^^'^ is the Dp-brane tension, are worldvolume coordinates 
and gab is the induced metric on the Dp-brane, Eq. ( |2.5D . The Dp-brane metric depends only 



on p (or z) and hence det (gab) will depend only on p (or z), so we will rescale the Dp-brane 



■^The full Dp-brane action of course describes the dynamics of a (7(1) worldvolume gauge field and all the 
worldvolume scalars (including 9{p), but also any others), and includes Wess-Zumino terms that describe the 
coupling of the Dp-brane fields to background fields, namely Ramond-Ramond form fields and the Neveu- 
Schwarz B-field. The equations of motion allow us to set to zero all components of the worldvolume field 
strength, as well as any other worldvolume scalars (besides our 0{p)), in which case we may safely ignore the 



Wess-Zumino couplings and work with the action in Eq. (2.6) 
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action and work with an action density. Among the coordinates one wih be the radial 
direction. Another wih be the time direction. Integration over this direction simply produces 
a factor of 1/T. The Dp-brane additionally has the internal directions of the C S^. 
Integration over these directions will produce the volume of a unit-radius j-sphere, Vj. For 
example, the D7-brane has V3 = 27r^. The remaining directions are spatial directions inside 
AdSi- In 5^ X 5'^ slicing, these are the directions of the 5*"^ C S^. Integration over these 
directions then produces a factor of Vi-2- In the various other slicings, these directions are 
not all compact, so integration over these directions will not always yield a finite volume. 
We will divide both sides of Eq. ( |2.6D by the factor of 1 /T and the integral over the spatial 
directions inside AdSi, and define an action density 



To explain the boundary conditions on the worldvolume scalar, we begin in the high- 
temperature phase, where the background geometry is global ^(iS'5-Schwarzschild times S^. 
We begin here because the story is essentially the same as for Dp-branes in Poincare-patch 
74(iS'5-Schwarzschild times S^. In the high-temperature phase, two classes of Dp-brane embed- 
ding are possible. In the first class, which we will call "Minkowski" embeddings, the C 
that the Dp-brane wraps is, at the boundary of AdS'5-Schwarzschild, the maximum-volume 
equatorial . As the Dp-brane extends into jddSs-Schwarzschild, the volume of the shrinks 
(as described by 9{z)) and, indeed, may shrink to zero volume at some radius outside the 
horizon, z = z < zh ■ The Dp-brane does not extend past z, rather, it appears to end at z [6] 
and does not intersect the ^dS's-Schwarzschild horizon. Explicitly, the boundary conditions 
on 6{z) are that 9{z) = ^, and 6'{z) = cxd in order to avoid a conical singularity [13]. In the 
second class of embeddings, which we will call "black hole" embeddings, the C shrinks 
as one moves into the bulk of ^dS's-Schwarzschild but never collapses to zero volume, so the 
Dp-brane intersects the horizon, thus developing a worldvolume horizon. Explicitly, we have 
9{zh) € [0, ^), and 0'{zh) = in order for the embedding to be static. These two classes of 
embedding are distinguished by topology: in Minkowski embeddings the collapses to zero 
volume, while in black hole embeddings only the thermal circle collapses to zero volume, at 
the horizon. The "critical solution" has z = zh, that is, the critical Dp-brane ends precisely 
at the horizon, and both the and the thermal circle collapse to zero volume. 

As first explained in Ref. [13], analogous embeddings exist in the low-temperature phase, 
but now the center of thermal AdS plays the role of the horizon. Minkowski embeddings 
now become "branes ending away from the center," that is, the C collapses to zero 
volume for some nonzero p, with boundary conditions 9{p) = ^ and 9'{p) = 00. Black hole 
embeddings become "branes that reach the center," that is, the C never collapses 
to zero volume, and the Dp-brane extends all the way to /? = 0, with boundary conditions 
9{0) G [0, ^) and ^'(0) = 0. These two classes of embeddings are again distinguished by 
topology: for branes ending away from the center, the C collapses to zero volume, 




(2.7) 



where Mop = ^f'^Dpyj- From now on, we will refer to Sop as the Dp-brane action. 
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Figure 1: The two topologically distinct Dp-brane embeddings in thermal AdS. The vertical circle 
represents the factor inside the S^, which can shrink to zero size in the left configuration, the 
horizontal circle represents the S*'"^ factor inside of AdSi, which can shrink to zero size in the right 
configuration. 



while for branes that reach the center, the 5"*"^ C AdSi collapses to zero volume. These 
embeddings are depicted schematically in figure ^. The critical solution is now a Dp-brane 
that ends precisely at the center, so 6{0) = ^ and both the and collapse to zero 
volume at p = 0. 

We will now explain in detail how the mass m of the flavor fields in the SYM theory, and 
the thermal expectation value of (the super symmetric completion of) their mass operator, 
may be extracted from the Dp-brane geometry. To be concrete, consider the low-temperature 
phase. Inserting the worldvolume metric Eq. ( |2.5D into the action Eq. (|2.7l ), we find 



Sdp = J^DpJ dp coshp (sinhp)^-2(cos6')Vl + ^'2- (2.8) 

If we expand this action in a Taylor series in 0{p), we may identify, from the 0{p)'^ term, 
a mass-squared for the worldvolume scalar, M| = —j (in units of the AdS radius). Probe 
Dp-branes fall into two classes depending on whether or not M| saturates the Breitenlohner- 
Preedman bound [23], > -i(i-l)2, for a field in AdSi with mass-squared M^. The 
field 9{p) for probe D5-branes and D7-branes does not saturate the Breitenlohner-Freedman 
bound, while 0{p) for probe D3-branes does saturate the Breitenlohner-Freedman bound. 

One may relate to the dimension A of the SYM theory operator dual to 9{p) in the 
usual fashion, M| = [A — (i— 1)]A. What is the dual operator? The operator dual to 9{p) is 
given by the variation with respect to m of the SYM theory Lagrangian. We will denote this 
operator as Om- For example, consider the probe D7-brane. The dual field theory in this case 
is = 4 SYM theory coupled to M = 2 supersymmetric hypermultiplets in the fundamental 
representation of the gauge group. The operator Om is the sum of three terms: the mass 
operator for the hypermultiplet fermions, m times the mass operator of the hypermultiplet 
scalars, and a term coupling these scalars to an adjoint scalar. The exact operator is written 
in Ref. [24]. The subscript on Om is a reminder that this operator depends explicitly on m. 
For our purposes thinking of Om as the mass operator of flavor fermions will be sufficient (we 
will not need the explicit form). 



The DBI action, Eq. (2.8), yields the equation of motion for 9{p). Given a solution, we can 
extract the value of m from the asymptotic behavior of 9{p) [2,3]. The asymptotic behavior is 
most easily studied (in both the low- and high-temperature phases) using Fefferman-Graham 
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coordinates, 

e{z) = z'-^-^ {^(0) + z 0(1) + • • • + [e(^2A~i+i) + ^{2A-*+i) log z] + • • • } . (2.9) 

For Dp-branes that do not saturate the Breitenlohner-Freedman bound, the leading term, 
with coefficient 0(o)! is non-normalizable, while the sub-leading term with coefficient 0(2A-i+i) 
is normalizable. The equation of motion for 6(z) determines all other coefficients in terms 
of these two. These two coefficients will be fixed by the boundary conditions explained 
above, which thus completely specify the asymptotic behavior. In particular, for these cases 
the coefficient 0(q) is related to the mass of the fundamental-representation fields as m = 
0(o)/(27ra'), or equivalently 6'(o) = 2Tx\f\m. Notice that 6'(o) is dimensionless and that the 
right-hand side of this equation is written in units of the AdS radius. From the SYM theory 
perspective, in x slicing, the more natural scale is the radius of the 3-sphere, i?3. In 
our units, i?3 = 1/2, so we may write 6(o) = 47rmi?3/\/A. 

For Dp-branes that do saturate the Breitenlohner-Freedman bound, which have 2A = 
i— 1, the coefficients 0(o) ^^^d 0(2A-j+i) ^ire identical. The leading logarithmic term, with 
coefficient ''A(2A-i+i) = V'(o)' is now non-normalizable, while the term with coefficient 0(o) 
is normalizable. In the general analysis of section ^, we will use the notation appropriate 
for Dp-branes that do not saturate the Breitenlohner-Freedman bound, namely 0(o) as the 
coefficient of the non-normalizable term and 0(2A-i+i) ^-s the coefficient of the normalizable 
term. Converting to Dp-branes that do saturate the Breitenloner-Freedman bound is easy: 
simply replace 6'(o) V(o) and i9(2A-i+i) ^ 6'(o) • 

Following Ref. [25], we can extract the thermal expectation value (Cm) from the asymp- 
totic coefficients using holographic renormalization. Let F denote the contribution that the 
flavor fields make to the free energy density. (Cm) is given by II . Using the AdS/CFT 
correspondence, F is given by the value of the on-shell DBI action Eq. (|2.7D , so we need to 
take the variation of the on-shell DBI action with respect to the leading, non-normalizable 
asymptotic value of Oi^z). The problem that arises in doing this is that the radial integration 
will diverge at the boundary, z = 0. In holographic renormalization, we introduce a regulator 
by cutting off the integration at some small value, i.e., we integrate only to z = e. We then 
introduce counterterms on the z = e hypersurface to cancel divergences before removing the 
regulator by sending e — > 0. These counterterms are written in terms of 0(e) and the induced 
metric on the z = e hypersurface, whose determinant we denote as 7. We denote the regulated 
action, plus counterterms, as Sreg-, so that S^eg = F ■ We then have, for Dp-branes that do 
not saturate the Breitenlohner-Freedman bound [25,26], 

For Dp-branes that do saturate the Breitenlohner-Freedman bound, we must make the re- 
placement e~^ e~^loge [25]. The counterterms and explicit formulae for (Om) for the 
D7-brane, D5-brane and D3-brane are given in Refs. [13,25]. 



- 10 - 



Let us illustrate the procedure with an example: the probe D7-brane, which is extended 
along AdS^ x 5^ and so has i = 5, j = 3 and M| = —j = —3. The D7-brane thus does not 
saturate the Breitenlohner-Freedman bound for AdS^, > —4. The dual operator Om has 
dimension A = 3. 9{z) has an asymptotic expansion 

e{z) = 0(0) Z + 0(2) + ^(2) zHogZ + -- - , (2.11) 

where ^^(2) is fixed in terms of 0(o) as '(/'(2) = l^^^(o)- The regulated action is Sreg = 
Sd7 + Y2k -^fc' where the counterterms are [13,25] 

Li = -l^, L2 = ^^n^, U = -U^oge)^(nlTV^ -\my (2.12a) 



4V M — 48V ' '^7' -^d — 32V^"6^/V / '^7 3 7^ 

U = \^e{t)\ U = -^^e{e)\ Le = ^log0{e)^n^e{ef. (2.12b) 

where for clarity we have suppressed a factor oi Mdi = NfT^fV^ = j^^XNjNc that appears 
in every counterterm. Here TZ]j and TZ^ are the Ricci tensor and Ricci scalar, respectively, 
of the induced metric on the z = e hypersurface. Using Eq. ( 2.10| ), we find that the thermal 
expectation value of Om is given by 

{Om) = Md7 [-20(2) + i^fo) + \^ ^(0) log ^(0)] • (2-13) 

We will also present numerical results for a probe D5-brane, which is extended along 
AdS^^ X S"^ and so has i = A and j = 2. The D5-brane has M| = —j = —2 and hence does 
not saturate the Breitenlohner-Freedman bound for — — f • The dual operator Om 

has dimension A = 2. In this case, the coefficient of the sub-leading, normalizable term is 
^(2A-i+i) = ^(1) • We will not present the counterterms explicitly. The result for the thermal 
expectation value of Om is (Om) = —J^DbO{i) [25], where Mob = NfT£)ry2 = ^^/X NjNc- 

3. Scaling Analysis 

In this section we perform a scaling analysis similar to that of Refs. [11, 14, 27] for probe 
Dp-branes in global thermal AdS, probe Dp-branes in global ^ddS'-Schwarzschild, and probe 
D7-branes in the various other slicings of global thermal AdS. We find the relevant scaling 
exponents for probe Dp-brane embeddings and determine when those exponents are real or 
complex. We also confirm via numerical analysis that complex exponents signal a first order 
transition in the dual SYM theory. 

3.1 Probe Dp-branes in x S"^ Slicing: Low-temperature Phase 

We begin with probe Dp-branes in global thermal AdS, corresponding to the SYM theory in 
the low-temperature, confined phase. We first note that the flavor physics, within the low 
temperature phase of the SYM theory, is completely independent of temperature. The global 



thermal AdS metric Eq. (2.1) does not depend on the temperature T, which thus appears 



in the DEI action Eq. (2^) only via an overall factor from integration over the thermal S^ 
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(which is then hidden in the rescaled action Eq. ( |2.7D ). Consequently, the equation of motion 
for 0{p), and its boundary conditions, are T-independent. As a result, the expectation value 
{Om) in the dual SYM theory will be independent of T. This temperature independence is 
not unexpected — it is a property of the leading large- A'^c behavior of generic observables (not 
involving Wilson loops which wrap the thermal circle) in the low temperature phase of non- 
Abelian gauge theories [28-30]. Our supergravity analysis is thus valid for any temperature 
below the Hawking-Page transition temperature, T < Thp = The remaining relevant 
scales in the SYM theory are the hypermultiplet mass m and the radius Rj, of the 3-sphere. 
The physics we are interested in will depend only on these scales, or more accurately on their 
dimensionless product mR^. 

We now turn to the scaling analysis. We focus on the region near the center of AdS, and 
expand the metric to leading nontrivial order using sinhp ~ p, cosh /? ~ 1. We first consider 
the critical solution, that, is, fluctuations of the form 6(p) = tt/2 + 59{p) with S9{p) small. 
We may then use cos^ ~ —69. The critical solution has the boundary condition 69{0) = 0. 
The Dp-brane action, Eq. (|2.8D, becomes 



Sdp = Ndp j dp p'-^ 69^ Vl + S0'^ , (3.1) 
with, once again, i + j = p + 1. The resulting equation of motion for 66 (p) is 

p6e69" + [{i-2)69 6e' -jp] {1 + 66'^) = 0. (3.2) 
This equation has an important scaling symmetry under which 

69{p) p69{p) and p ^ pp , (3.3) 

for real, positive p. In other words, a single solution 69{p) = f{p) gives rise to a one-parameter 
family of solutions 69 (p) = p^^ f{pp). The solution of Eq. ( |3.2| ) for the critical embedding, 
which we denote 69* (p), is 

^(^*ip)^P\f^- (3-4) 

Notice that 69* (p) is invariant under the scaling transformation. 

Next, we seek solutions near the critical one, still in the near-center region, of the form 

69{p) = 69*{p)+ap). (3.5) 
The fluctuation ^(p) has the linearized equation of motion 

p^^" + ap^' + a^ = 0, (3.6) 
with a = j + {i — 2). This has solutions ^{p) = ibp^± with 



/3± = i (1 - a) ± -6a + l . (3.7) 
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The (3± are the scahng exponents. Depending on the value of a, the scahng exponents P± 
may be real or complex. Physically, a = j + {i — 2) is the dimension of the sub-manifold of the 
Dp-brane worldvolume that collapses to zero volume in the critical solution: the j dimensions 
of the C and the i—2 dimensions of the C AdSi. For the present case of probe 
Dp-branes in x slicing, the dimension a of the collapsing submanifold is simply p—1, 
but this will not be the case in subsequent examples. 

Explicit values of the exponents, for various choices of a, ap- 
pear in Table |^. As shown in the table, the critical value of a is 
6: for a < 6 the near-center, near-critical solutions have complex 
exponents, while for a > 6 they have real exponents. 

The perturbed solution, 



5e{p) = 6e*{p) + a. 



(3.8) 



a 


/3± 


7 


-3± ^2 


6 


_5 ± 1 


5 


-2±i 


4 


2 ^ ^ 2 


3 


-l±iV2 


2 


2 ^ ^ 2 



must obey the scaling symmetry. This implies that the coefficients 
a± must scale as 

a± /i^-^± a± . (3.9) 

The coefficients a± are determined by the boundary conditions, 
specifically, the value of p or ^(0). This is why a± are "charged" 
under the scaling transformation: because p and 9(0) are not in- 
variant under scaling. 

A solution that is near the critical solution in the near-center 
region will remain so all the way out to the asymptotic region.^ 
The boundary condition deep inside AdSi (the value of p or 6^(0)) 
must fix both the asymptotic coefficients 0(o) and 9(2A-i+i) and the 

near-center coefficients a± . The asymptotic coefficients can thus be thought of as functions of 
the near-center coefficients. Sufficiently close to the critical solution, the a± will be very small 
and hence the asymptotic coefficients may be linearly related to the near-center coefficients. 
Put more simply, we may Taylor expand 9(^q^ and ^{2A-i+i) to linear order in a±. Let ^^q^ 
and 9'^2A-i+i) denote the asymptotic coefficients of the critical solution. We then have, for a 
near-critical solution, 



Table 1: Scaling ex- 
ponents P± for probe Dp- 
brancs whose critical solu- 
tion has a collapsing sub- 
manifold of dimension a. 



^(0) - ^(0) = ^(0) «+ + ^(0) «- ' (3.10a) 
&{2A-i+l) - ^*2A-i+l) = ^(2A-i+l) + ^{2A-i+l) ' (3.10b) 

for some set of coefficients A^. Notice that the overall sign of ^{p) is arbitrary, so the form of 

/j-Qj aiiia i/^oj ^ yy^Qy that Is, for Dp-braues that end at the 



ho) < %) and 6'(o) > 



Eq. ( 3.10| ) is valid for both 
center and for Dp-branes that end away from the center. We may now ask what the scaling 
transformation of the a±, Eq. (^.9|), teaches us about the asymptotic coefficients, and what 
we may then conclude in the SYM theory about the behavior of {Om) as a function of m. 

^More formally, we expect the critical solution to be an attractor solution as one moves toward the AdS 
boundary. In Ref. [27], this was shown explicitly for Dp-branes in Poincare-patch yldS'-Schwarzschild. 
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For complex P±, if we perform the scaling transformation Eq. ( p.9| ) we find that, as func- 
tions of n, (6'(o) — ^(0)) ^^'^ (^(2A-j+i) ~ ^(2A-i+i)^ ^^^^ acquire terms of the form ^i-Re/3± 
times either a sine or cosine of Im.p± log^. Notice that 1 — Jie l3± > in all cases. Conse- 
quently, sending /i ^ produces solutions which, in the (^(o); ^(2A-i+i)) plane, spiral inward 
toward a limit point which corresponds to the critical solution. More precisely, the arbitrary 
overall sign of (,{p) leads to two intertwined spirals, as illustrated below in Fig. ^(b). 

In the dual field theory, we may argue, following Refs. [11, 14], that {Om) must jump 
discontinuously from the outer branch of one spiral to the outer branch of the other as the 
mass m is varied. Such discontinuous behavior signals a first-order phase transition. We 
then expect discontinuous behavior in generic observables associated with the flavor fields, of 
which {Om) is simply the most convenient to compute from the supergravity perspective. 

To be concrete, let us analyze the example of probe D5-branes extended along AdS/i x S'^, 
for which a = 4 and P± = — | ± i^. Recall that in this case m = 6'(o)/(27ra'), and (Om) = 
— 7Vd50(i) where the operator Om has dimension A = 2 and Md^ = '^^^^'^^f^c- Figure ^ 
shows the value of Sreg/-^D5 for the Z)5-brane as a function of 0(o), and figure ^ shows the 
numerical result for solutions in the (6'{o)5^(i)) plane. 




(a) (b) 

Figure 2: (a.) Sreg/J^D5 as a function of 0(0) for the D5-brane probe in global thermal AdS in 
X slicing. The red curves correspond to D5-branes that reach the center of AdS while the black 
curves correspond to D5-branes that end away from the center. The point at which the red and black 
curves meet corresponds to the critical solution, (b.) Close-up of (a.) near the critical solution. The 
inset figure shows a further close-up. The vertical line in the inset figure indicates where the transition 
occurs, at the critical value 6*^^^** « 2.402. 

In these and all subsequent plots for global thermal AdS, the red curves arise from Dp- 
brane solutions that reach the center of AdS, while the black curves arise from Dp-branes 
ending away from the center. For global Ad^s-Schwarzschild, red curves arise from Dp-branes 
that reach the horizon while black curves arise from Dp-branes that end outside the horizon. 
The critical solution will thus always be where the red and black curves meet. 

In figure we see that the free energy, which is always the minimal value of the regulated 
action Sreg, has a small kink (a discontinuous first derivative) as a function of 0(o)- Since 
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(a) (b) 

Figure 3: (a.) as a function of 6(Q-f for the D5-brane probe in global thermal AdS in x 
slicing. The red curves correspond to D5-branes that reach the center of AdS while the black curves 
correspond to D5-branes that end away from the center. The point at which the red and black curves 
meet corresponds to the critical solution, (b.) Close-up of (a.) near the critical solution. The vertical 
line indicates where the transition occurs, at the critical value S^q^'* w 2.400. As we increase 0(o), 
moving along the red curve from left to the right, the physical value of 9(i) = —{Om)/J^D5 jumps 
upward at 0?^*, from the red curve to the top-most arm of the black curve. 



m = ^(0) / (2'/ra'), we conclude that a first-order phase transition occurs in the SYM theory as a 
function of the mass m. (The free energy is always continuous, but at a first order transition its 
first derivative jumps.) The transition occurs at the critical value^ = ATimR^/y/X = 2.400 
or equivalently m = VXOf^f /{AttRs) = 0.191 VX/R^. Using (Om) = |^ = -J^DbO^i), we see 
precisely this behavior in figure as jumps discontinuously from one arm of the spiral 
to the other. When the embedding equations have multiple solutions for a given value of 
0{o), only those solutions that minimize the free energy (Sreg) represent genuine equilibrium 
states. These minimal free-energy solutions only lie on the outermost branches of the two 
spiral arms. As shown in Refs. [14, 19], moving inward along the red or black curves toward 
the critical solution, a new tachyon appears in the meson spectrum at every turn of the spiral, 
providing a clear signal of instability. Notice also that m* ~ I/-R3, as expected from the fact 
that generic observables will be T-independent in the low-temperature phase, as mentioned 
above. Clearly this transition is a finite-volume effect. 

When the scaling exponents f3± are real, the story is very different. Now the scaling trans- 
formation of Eq. ( ^ ) is a simple rescaling, and no spiral will appear in the (^(o); ^(2A-i+i)) 
plane. The most interesting example is the probe D7-brane,^ extended along AdS^ x S^, 

^Notice that the result for Sjoj* as computed from the free energy in figure ^ is ~ 2.402, which is 

slightly larger than the value ~ 2.400 obtained from figure where we used an equal-area method. We 
attribute the difference to numerical error in estimating the location of the "kink" in figure ^. Here, and for 
all subsequent branes, the location of the transition that we will use is that computed from the equal-area 
method. 

^We have also formally analyzed a D8-brane extended along AdSs x S"*, which has a = 7, and found 
behavior similar to that for the D7-brane. A D8-brane carries no conserved charge in type IIB supergravity, 
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(a) (b) 

Figure 4: (a.) Sregl-N'm as a function of 6'(o) for the D7-brane probe in global thermal AdS in S*^ x 
slicing, (b.) Close-up of (a.) near the critical solution. The critical solution (where the red and black 
curves meet) has 91^^-, ~ 2.198. 




(a) (b) 

Figure 5: (a.) 6(^2) as a function of O^q-^ for the D7-brane probe in global thermal AdS in x 
slicing, (b.) Close-up of (a.) near the critical solution. The critical solution (where the red and black 
curves meet) has 91^-, w 2.198 and 91^^ ~ 5.230. 



for which a = 6 and f3± = — | it ^ = {—2, —3}. Recall that in this case ra = 6 (q-^/ [211 a') 
and (O m) is given by Eq. ( p.l3D where the operator Om has dimension A — 3. Recall also 



that Nd7 = (2^)4 XNfNc- Figure ^ shows Sreg/J^m as a function of 0(o) the D7-brane, 
which has a continuous first derivative. Figure ^ shows 6(2) as a function of 0(o)) which is 
single- valued. The critical solution has 0^q-j ~ 2.198 and 0*2) ~ 5.230. In terms of SYM the- 
ory quantities, the critical solution has m* = \/A 0*Q)/(47ri?3) = 0.1748\/A/i?3. (Notice again 
that m* ~ l/iis.) We conclude that no first-order transition occurs in the SYM theory. We 
expect some kind of non-analyticity in the SYM theory, however, because in the supergravity 
description the topology of the D7-brane changes in passing through the critical solution. In 
section Q we will show that certain observables have divergent behavior at the critical point, 
and we will examine how this transition affects the meson spectrum. 



however, and is therefore unstable. 
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3.2 Probe Dp-branes in x Slicing: High-temperature Phase 

We will now perform the scaling analysis for probe Dp-branes in global yldSs-Schwarzschild, 
corresponding to the SYM theory in the high-temperature, deconfined phase. Our result 
will be essentially the same as that of Refs. [11, 14]: all probe Dp-branes have complex 
exponents P± and hence exhibit a spiral in the (^(o)i ^(2A-i+i)) plane, implying that first- 
order transitions appear in the respective SYM theories. For the D7-brane, this must be 
the case in order to agree with the numerical results in x slicing [13] and in Poincare 
patch AdSs-Schwarzschild [7-9,11-14], which we interpret as the high-T limit of the theory 



formulated on global ^(iS'5-Schwarzschild, as explained in section 2.1. 

The critical solution is now a Dp-brane that ends precisely at the horizon. At the horizon 
only the thermal inside AdS'j-Schwarzschild collapses to zero volume. This is the key 
difference from thermal AdSi. For the critical solution, the C will still collapse, but 
now all probe Dp-brane critical embeddings will have the same collapsing S^. The 5*-^ that 
collapsed in thermal AdSi are replaced with this in AdS'j-Schwarzschild, so we expect 
a = j + 1. We can thus jump to the answer: as the largest value of j that allows for 9{z) is 
j = 4, none of the Dp-brane probes can have a > 6, and hence all must have complex (3±. 

To confirm this, we proceed in the same spirit as above. We now focus on the near- 
horizon region. Let z = zh + Z and 9{z) = f + S9{Z). We expand the metric coefficients to 
leading nontrivial order. The induced Dp-brane metric is then 

2 

dslp = (1 + 69'^) dz^ + — — dT^ + :f{zh) dnf_2 + <5^^ dnl (3-11) 

The Dp-brane action, ignoring overall Z-independent constants that do not affect 56^s equa- 
tion of motion, is 



Sopoc j dZZ 5eW^ + 50'^ . (3.12) 



This is of precisely the same form as Eq. ( p.l[) but with i—2 1 so indeed a = j+\ and all T>p- 
brane probes will have complex exponents j3± and exhibit a spiral in the (^(o)) ^(2A-j+i)) plane. 
The action Eq. is in fact identical to the near- horizon action written in Refs. [11, 14] 



for Dp-brane probes in the Rindler space that arises as the near-horizon geometry of Poincare 
patch AdS'-Schwarzschild. 

3.3 Probe D7-brane in Other Slicings: Low- Temperature Phase 

By using different slicings, leading to different boundary geometries, we can change the di- 
mension of the critical solution's collapsing submanifold, as different slicings lead to spheres of 
different dimension, 5' C AdS^ for I = 0, . . . , 4, that collapse to zero volume at the center of 
AdS^ (see Eq. (|2.3D). We will focus on the D7-brane and the low-temperature phase because 
lower-dimensional Dp-branes, or Dp-branes in the high-temperature phase, will not be able to 
reach a > 6. For the D7-brane, the dual SYM theory is AA = 4 SYM theory coupled to mas- 
sive M = 2 hypermultiplets, formulated on different four-manifolds, in the low-temperature, 
confining phase. 
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Using the AdS^^i x slicing of AdS^, the induced D7-brane metric is 
dsl)j = [1 + 9'{pf] dp^ + cosh^p ds\j^s^_^ + sinh^/? dVLf + cos^ e{p) dflj , 
and the D7-brane action is 

Soj=Md7 j dp (coshp)^-' (sinhp)'(cos0)Vl + ^'2. 
In the near-center hmit this becomes 

Sjj^=Mm I dpp^ 56'^ ^/ 1 + 56'^ 



(3.13) 



(3.14) 



(3.15) 



which is the same as Eq. ( p.l[ ), but with i — 2 — > / and j = 3, so a = Z + 3 and we will have 
complex exponents I3± for Z = 0, 1, 2 and real exponents j3± for / = 3, 4. 

As an example of complex exponents, consider the / = 2 case, AdS2 x 5^ slicing, with 
a = 5 and complex exponents /3± = — 2ibz. We expect a spiral in the (^(o))^(2)) plane. Figure 
P shows the numerical result for 9^2) as a function of 0(o)) which indeed exhibits a spiral, so 
we again have a first-order transition. In this case, the transition occurs at the critical value 



m/(27ra') = 1.6557 or equivalently m = y/XOpI^ / {2tt) = 0.2635 \/A times the inverse 



ncrit 
^(0) 

of the AdS curvature radius. 




(a) 



(b) 



Figure 6: (a.) 9(^2) as a function of ^(q) for the D7-brane in AdS2 x S'^-sliced thermal AdS^. (b.) 
Close-up of (a.) near the critical solution. The vertical line indicates where the transition occurs, at 
the critical value 6*^^]* « 1.6557. 



We have two examples of real exponents. The first example is the I = 3 case, x 5^ 
slicing, which was examined above (see figures ^ and H) and will be examined in detail in 
the next section. The second example is the I = 4 case, slicing, with a = 7, giving real 
exponents, (5± = —3 ± \/2. We expect 9^2) to be single- valued as a function of ^(q). Figure 
^ shows the numerical result for 9(^2) as a function of 0(o)i which is indeed single- valued. 
We have verified numerically that the solutions exhibit scaling behavior with the exponents 
P± = -3± V2. In this case the critical solution has 9*^^ ~ 2.6395 and 0^2) ~ 14.5943. In 
terms of SYM theory quantities the critical solution has m* = \/X0*Q^/(27r) ~ 0.42-v/A times 
the inverse of the AdS curvature radius. 
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(a) (b) 

Figure 7: (a.) 9(2) as a function of O^q^ for the D7-brane in S"*-sliced thermal AdS. (b.) Close-up 
of (a.) near the critical solution. The critical solution (where the red and black curves meet) has 
er„. w 2.6395 and 6*. « 14.5943. 



4. Characterizing the Transition 

Let us summarize our story so far. Whether a Dp-brane probe exhibits a spiral in the 
(^(0)' ^(2A-i-i-i)) plane or not depends upon the value of a, the dimension of the submanifold 
of the critical embedding that collapses to zero volume. When a < 6, a spiral does appear 
and, as discussed above, we may conclude that the dual SYM theory has a first order phase 
transition associated with the flavor fields. 

To study what happens when a > 6, we will restrict our attention to the D7-brane in 
X slicing. The D7-brane undergoes a topology change, but 0(2) is single-valued as a 
function of 0(o) (o^ mass). In the SYM theory, we expect a transition to occur, but evidently 
it is not first order. 

In this section, we will exhibit divergences at the critical point in several observables in 
the dual SYM theory, and examine the associated critical exponents. We will also show that 
the transition manifests itself in the meson spectrum as a cusp in the spectrum of a particular 
scalar meson. We will also argue that the transition will be "smoothed out" for any large but 
finite value of A, or in other words that the transition is an artifact of the A — > oo limit. 

4.1 The Behavior of (O^) 

We want to know what the scaling transformation Eq. ( p. 9] ) implies for the observable (Om)- 
Returning to Eq. (3.10), we perform the scaling transformation in Eq. ( |3.9| ), using 1 — /3+ =3 



and 1 — /?_ =4, with the result 



^(0) - = ^(0) "+ /^'^ + ^(0) ' (^-1^) 

0(2) - 0*2) = ^ J) "+ + ^(2) • 

Notice again that, as in Eq. ( |3.10| ), the form of these equations (and the equations below that 
follow from them) is the same for both 0(o) > ^fo) ^(0) ^ ^fo)' Using the first equation 
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to eliminate fx gives, to leading order in the deviation from criticality, fi ~ (0(o) ~ ^('o))^^^' 
Plugging this into the second equation will generate a non-analytic (^(q) — ^('o))^''^ term as 
well as an analytic (6'(o) — ^?q-,) piece. Translating to field theory quantities using 0(o) 



oc m 



and Eq. ( 2.13| ), we find that (Cm) contains non- analytic terms of the form 

(O™) ~(m-m*)^/3^--- , (4.2) 

in addition to contributions analytic in m near m* . Here • • • is an expansion in higher powers 
of (m — m*)^/^. The key point is that (Om) has a divergent second derivative with respect 
to m at m* . The expectation (Om) is the derivative of the free energy F with respect to m, 
and hence the third derivative of the free energy diverges as m — > m* , 

F~ (m-m*)-2/3, (4.3) 



from which we identify a critical exponent of 2/3. 
4.2 Static Test Charges 

Consider, in the dual SYM theory, the interaction between static fundamental-representation 
test charges. Focus, for simplicity, on the case of an infinitely heavy quark and antiquark which 
are maximally separated — sitting at opposite poles of the S^. The change in free energy 
due to inserting the static quark and antiquark is given by {—T times) the logarithm of the 
expectation value of the Polyakov loop (or Wilson line) correlator, with the two Polyakov 
loops sitting at antipodal poles in space. 

Roughly speaking, in the dual gravitational description (in the Nc ^ oo and A ^ oo 
limits), this Polyakov loop correlator is given by the regularized minimal area of a string 
worldsheet with an x boundary, where the 5^'s wrap the thermal circle and are max- 
imally separated in the boundary spacetime. More precisely, the objects that have simple 
supergravity descriptions are the supersymmetric extensions of Wilson loops ( "Maldacena 
loops"). These correspond to adding a line integral of a linear combination of the SYM 
scalars to the exponent of the Polyakov loop [31] . Exactly what linear combination is deter- 
mined by choosing, arbitrarily, some six-dimensional unit vector, or equivalently some point 
in S^. In the dual description, the boundary of the string worldsheet must be located at the 
chosen point in the factor of the boundary geometry. 

We will focus, for simplicity, on the maximally symmetric case where the Maldacena loops 
are located at a pole of the S^, or in other words are maximally separated (in the boundary 
geometry) from the D7-brane wrapping an equator of the S^. This choice preserves the 
full SO{4) /^-symmetry of the SYM theory with fundamental hypermultiplets. The antipodal 
locations of the loops means that the correlator is also invariant under an SO (3) subgroup of 
the 50(4) spatial rotation symmetry group, as well as U{1) time translations. 

Given this setup, an extremal string worldsheet, preserving all the symmetries of the 
correlator, will be one that stretches across global AdS^ while everywhere sitting at the pole 
of the S^. This is illustrated in figure ^ (B.). 
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Figure 8: (A.) The disk depicts global AdS^ in x slicing. The perimeter represents the boundary 
while the center represents the center of AdS^. The shaded/striped area represents a D7-brane that 
ends before reaching the center. The red lines represent two segments of a string, each stretching from 
the boundary to the D7-brane's endpoint. Each of the three circles between (A.) and (B.) represents 
the S^. The D7-brane, wrapping an C S^, is depicted as a black line. The location of the string(s) 
is depicted by a red dot. At the boundary, the D7-brane wraps the equatorial S^, but at its endpoint 
it reaches the pole of the S^, where the strings also sits (of the circles between (A.) and (B.), compare 
the top one and the lower left one). (B.) The same picture as in (A.), but now the D7-brane does not 
end at finite radial coordinate, rather it extends all the way to the center. Indeed, the figure depicts 
the trivial embedding, in which the D7-brane wraps the equatorial everywhere. Notice that in this 
case the string and D7-brane never coincide on the S^. 

If the D7-brane embedding fills the the whole AdS^ space and reaches the center, then 
the wrapped C never collapses to zero volume (so < ^(p = 0) < and the D7-brane 
never intersects the string worldsheet. This is not the case for D7-brane embeddings that end 
outside the center (at some p = p where 6{p) = ^). Such D7-branes reach the pole of the 
internal and intersect the string worldsheet at p. This allows the string to break, so that 
its worldsheet ends on the D7-brane instead of extending all the way to the center of global 
AdS. This is illustrated in figure ^ (A.). Phrased differently, for Minkowski embeddings of 
the D7-brane, a maximal-symmetry string worldsheet with an (Si)^ boundary exists — with 
two of the S^^s on the AdS boundary and the other two on the D7-brane at a point of closest 
approach to the global AdS center. Again, the two possibilities are depicted in figure ^ 

The potential energy of the static quark and antiquark is directly proportional to the 
length of the string. The difference in length of the strings in the two configurations is just 
the diameter of the "hole" in the center of the D7-brane configuration. For the near-critical 
embeddings, this is directly proportional to the scaling parameter /i. Therefore the difference 
in the static quark-antiquark potential energy (for this maximally-symmetric configuration) 
between near-critical and critical values of the quark mass will proportional to /i ~ (m— m*)^/^ 
(or (m* — m)^/'^ for m < m*) — and thus non-analytic in the quark mass at m* . 

The behavior of this Polyakov loop correlator has an intuitive explanation in terms of 
screening lengths. To understand this, first consider the physics of QCD (in flat space) at low 
temperature, where two infinitely massive test quarks will have a fiux tube, or QCD string, 
stretching between them. The energy E in this QCD string is its tension, o", times its length. 
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L: E = aL. The theory also has dynamical quarks, however, which can cause sufficiently 
long flux tubes to break. If the test quarks are pulled sufficiently far apart, then the energy 
in the string will exceed the point at which the creation of a dynamical quark-antiquark pair 
is energetically favorable. In other words, when E > 2m for a dynamical quark mass m, the 
string can snap and the single heavy-heavy meson will decay into two heavy-light mesons. 
Define a "screening length" as the length Lg at which the pair production occurs and the 
QCD string snaps. Equating aLs = 2m, we see that the screening length is proportional to 
the mass of the quark, Lg oc m. 

Now consider the M = 4, SYM theory formulated on coupled to massive M = 2 hy- 
per multiplets. The story is a bit different since the theory is conformal. This implies that 
the static quark potential is Coulombic, and vanishes at large separation. Naively, one might 
think that this would imply that pair production of dynamical massive quarks would never 
be energetically favorable. (After all, hydrogen is not unstable to electron/positron pair cre- 
ation!) This argument relies on our weafc-coupling intuition, however. If the coupling is 
sufficiently strong, then the light-heavy binding energy can become comparable to the light 
quark mass m. In this regime a heavy-heavy meson, pulled apart to large separation, can 
indeed decay to a pair of light-heavy mesons via pair production of a light quark-antiquark 
pair. This scenario is sometimes called Gribov confinement [32]. For the AA = 4 SYM theory, 
this process was analyzed in detail in Ref. [33]. The screening length, beyond which the 
system will pair produce dynamical quarks, must by scale invariance be proportional to the 
inverse of the light quark mass, Lg oc 1/m. 

Finally, consider the M = 4 SYM theory formulated on with M =2 matter. As we lower 
the dynamical quark mass we find that at some point this screening length becomes larger 
than the diameter of the sphere, Lg > 2R^. Beyond this point, a (maximally separated) static 
quark-antiquark pair is stable against pair production of dynamical quarks. Our Polyakov 
loop correlator exhibits precisely this behavior. 

4.3 The Meson Spectrum 

In the = 4 SYM theory coupled to A/" = 2 matter on x M'^, the first-order phase transition 
of the fundamental-representation fields is characterized by ''meson melting."' At low tem- 
perature the meson spectrum is gapped and discrete while at high temperature it becomes 
gapless and continuous [11, 14,34]. A natural question for the theory on S^, then, is what 
happens to the meson spectrum at the transition we have found? 

For conceptual simplicity, and for technical reasons that will become clear shortly, in 
this section only we will take the temperature T to be precisely zero. We are allowed to 
do so because, as explained above, the D7-brane physics we are interested in is independent 
of temperature in the low-temperature phase. We are thus free to take T = and study 
the SYM theory formulated on the four-manifold M x 5^ with R the (Minkowski-signature, 
non-compact) time direction. 

At zero temperature we have a simple intuitive picture of what happens. Consider first 
the infinite- volume case at finite temperature. If the mass m is fixed and we heat the system 
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Figure 9: (a.) Scalar meson frequency squared (times tlie AdS radius squared) versus 6'(o). The 
critical solution has = 2.198 and (w*)^ — 5.7. (b.) Close-up of (a.) near the critical solution. Our 
analytic arguments (confirmed numerically in figure p^) show that the curves behave as {0(q) ~^(o))^^'^ 
near the critical solution. 



up, stable bound mesons fall apart at the "melting" temperature. What we are instead doing 
here may be viewed as fixing the mass and then squeezing the system, at zero temperature, 
into smaller and smaller volume and asking what happens to mesonic bound states as R3 
passes through the critical radius = 0.1748\/A/m. We expect the mesons to fall apart as 
their zero point energy due to confinement within the finite volume becomes larger than their 
binding energy. 

Before discussing the meson spectrum in detail, we should define precisely what we mean 
by a meson mass. On M x S^, a state is classified by its SO{4) angular momentum /(/ + 2), 
for non- negative integer /, and its energy eigenvalue w, governing the behavior under time 
translations. We will use this energy eigenvalue co to characterize a meson state, and will 
refer to to as the meson mass.^ We will only consider states that are in an Z = s-wave on S^, 
and compute the dependence of on the quark mass m, or more precisely (but equivalently) 
on 0(0)- States with higher / are expected to have higher energy. As explained in Ref. [35], 
the meson spectrum can be computed in the supergravity description by expanding the DBI 
action to quadratic order in fluctuations and then solving the resulting linearized equations of 
motion. The fluctuations may be those of either the embedding geometry or of the D7-brane 
worldvolume gauge fields. 

We first study a fluctuation of the D7-brane geometry of the form 

0{p)^0{p,t)=eBip) + <P{p)e'^', (4.4) 

®In some literature on holographic mesons in flat space, the name "meson mass" is used for the magnitude of 
an imaginary spatial wave vector k for which an eigenstate of the supergravity small fluctuation operator exists. 
These eigenvalues characterize the leading long distance faU-off of equilibrium Euclidean space correlators in the 
dual field theory. On 5*^, no (interesting) analogue of this definition exists. The spatial (angular) momentum 
simply takes values + 2) for non- negative integers I; since the space is finite we have no notion of the 
asymptotic behavior of correlation functions. 
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where Ob{p) is a static background solution, as computed in section 3.1. This fluctuation will 
be dual to a scalar meson in the SYM theory. To determine w, we employ a shooting technique 
in which we fix the boundary condition for 0(p) at the AdS boundary and iteratively adjust 
u) until the physical boundary condition at the D7-brane endpoint (or the center) is satisfied. 
At the AdS boundary, we demand that the fiuctuation be normalizable and hence must scale, 
in Fefferman-Graham coordinates, as 0{e^) with the cutoff z = e. For D7-branes that reach 
the center, we then demand that the fiuctuation must have vanishing first derivative at the 
center. For D7-branes ending away from the center, we require that the C collapse 



without a conical deficit, as explained in section Translating the results into field theory 
quantities, the resulting spectrum of versus 0(o) — 47r\/A Jn-Rs for this scalar meson is 
shown in figure We see a pronounced kink precisely at the critical solution, 0(o) = ^(o)- 
The kink does not extend down to uj"^ = 0; the minimal value is ~ 5.7 times the square of 
the AdS radius, or equivalently lo"^ 22.8i?|. 

Using scaling symmetry arguments, we can show analytically that this kink in the scalar 
meson spectrum has the precise form uj — uj* ^ (m — m*)^^^ as m ^ m* (from above or 
below). We begin by solving analytically for the fiuctuation in the near-center limit. At zero 
temperature the action for 0{p,t) is 

So7 = Mm j dp dt (cos ef (sinh pf\/{l + 0'2)(cosh2 p - 6^) + 9"^ 9'^ . (4.5) 



We take the same near-center limit as in section 3.1, insert Eq. ( |4.4j ), and expand to quadratic 
order in the fiuctuation (p{p) to find the linearized equation of motion. We choose the back- 
ground solution to be the critical solution 9b{p) = 9*(p) = ^ + p, so the frequency of the 
fiuctuation is the critical one, uj = uj*. The resulting equation for (/'(p), 

p^ f{p) + 6p (p'{p) + (6 + 2a;* V') <t>{p) = 0, (4.6) 



has a solution 



(t>{p) = \ 
P'^ 



* ( f5 * \ ^ * sin(\/2u;V) 
Cg cos(v2a; p) + Ci 



4p-' + cip-^ + oip-'), 



(4.7) 



where Cq and are integration constants. The background is the critical solution, for which 
the brane ends at the center, p = 0. Normalizability at p = requires Cq = 0. We could then 
fix the value of uj* by imposing normalizability at the AdS boundary as follows. 

Normalizability at the boundary requires that </'(p)'s contribution to the leading, non- 
normalizable asymptotic coefficient (which was ^(q) i'^ Fefferman-Graham coordinates) must 
vanish. We denote this contribution as 

<^*icl,u;*) = fi{u;*)cl + 0{cf) = 0, (4.8) 

where, as indicated, this may be a function of and uj* and we have linearized in c^, which 
is taken to be small. We could then solve the equation fi{uj*) = to find the value of u;* . 
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We next want to find the shifted frequency, uj = lo* + 5uj, of fluctuations of near-critical 
solutions. More precisely, we need to determine how 5uj = uj — oj* scales with ^ ~ (m — m*)^/'^. 
Near-critical solutions will have nonzero values for cq and ci = c\ + 5ci so that: 

e{p, t) = p + a^ + a+ p-^ + ^(p) e*"* (4.9) 
= p+{a-+co e*"*) + (a+ + ci e'"*) p~^ + 0{p-^) . 

We see that cq must scale the same way as a_, and ci the same as a+, that is cq ^ Ai'^co 
and ci — > /u^ ci. To fix the value of w, we again impose the condition of normalizability at the 
AdS boundary, which requires that the coefficient $(ci,co,u;) of the non-normalizable term 
vanish, 

= $(ci,co,u;) = /iMci+/2Mco (4.10) 
= [/i(.;*) + f[{u;*)6u;] (c^ + 5ci) + [^(w*) + f^{u*) 5u;] cq 
= fii^*) S^cl + 72(0;*) Co , 

where we have linearized everything treating 6ci, cq and 5uj as the same order of smallness 
and used /i(a;*) = 0. We may immediately solve for 5uj = to — uj* with the result 

* /2(w*) Co , . 

'^-'^ =-/[(:?) jj- 

Notice 72(0;*) and /{(w*) do not transform under the scaling symmetry: they are just numbers. 



Since co scales as p'^ while c^ scales as p^, the result of Eq. ( 4.11 ) shows that a; — ~ ~ 
(m — m*)^/^. Notice that these arguments hold for both m > m* or m < m* (the results do 
not depend on whether the near-critical solution in Eq. ( |4.9| ) has ^(o) > ^(o) ^(0) < ^(o))' 
Figure 10 shows our numerical data for In la; — w*! versus In |0(o) ~ ^ linear fit to the 

data, with good agreement between this asymptotic form and the numerical results. 

This non-analytic behavior appears to manifest itself only for the geometric fluctuation 
of the embedding function 9{p). Solving for the other fluctuation in the geometry (the other 
direction orthogonal to the S^) we can compute the spectrum for a second scalar meson. 
The result appears in figure |ll|. We see no kink at 9*^Qy We also computed the meson 
spectra corresponding to fluctuations of the D7-brane's worldvolume gauge field as described 



in Ref. [35]. These meson spectra are similar to figure 11. In particular, they exhibit no 
kinks. We confirmed that all of our meson spectra reduce to the known results of Ref. [35] in 
the limits of either zero or large mass. 

4.4 Light String States 



In section 4.1, we showed that the third derivative of the free energy F diverges at m* . 
In other words, some three-point coupling in an effective theory describing fiavored mesons 
in this strongly-coupled SYM theory is diverging as m approaches m* . One question this 
immediately raises is whether this power-law growth of the three-point coupling continues for 
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Figure 10: (a.) In |cj — w* | (with u in units of the inverse AdS radius) versus In [^(o) — 0*^^ | for D7- 
branes that end at the center, which have 0(o) < ^*o) (corresponding to the red curves in the previous 
figure). The solid black line is a numerical fit to a functional form Ci + C2 In |6'(o) — ^(0) I- O^^' analytic 
argument predicts C2 = 1/3. The numerical result is Ci ~ —0.38 and C2 ~ 0.35. (b.) The same 
quantities as in (a.) but now for D7-branes that end away from the center, which have 0(o) > (the 
black curves in the previous figure). The solid black line is a numerical fit of the same form as in (a.) 
with the result Ci « -0.24 and C2 ~ 0.36. 




Figure 11: Scalar meson frequency squared (times the AdS radius squared) versus 6'(o) for the meson 
holographically dual to the second fluctuation of the D7-brane geometry. The punchline here is that 
no "kink" appears, in contrast to the spectrum in figure ^ 

arbitrarily small values of m — m*, and if so, what is the physical meaning of this singularity? 
We will argue that for any large, but finite^ values of A, the scaling regime will be cut off at 
small values of m — m*, so the divergence in the three-point coupling is an artifact of the 
strict A ^ 00 limit. 

Our evidence comes from the string theory side of the correspondence: the scaling regime 
is cut off by stringy corrections. To see this, note that the scalar curvature of the D7-brane's 
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Figure 12: Log of minus tlie scalar curvature of the D7-brane, ln[— 7?._d7] (in units of the AdS 
radius), evaluated at the point of closest approach to the center of AdS, versus asymptotic coefficient 
6'(o) in 5'-^ X slicing. The red curve represents D7-branes that reach the center, for which ln[— 7?._d7] 
is evaluated at the center. The black curve represents D7-branes that end away from the center, for 
which ln[— 7?.£i7] is evaluated at the endpoint. 



induced metric, TZor, in the near-center limit, is 

66^69'^ 56 56' + 

For the critical embedding, 56* (p) = p, the curvature diverges as p~'^ near the endpoint 
p = 0. For near-critical solutions, the curvature is finite at the endpoint (or center), but 
grows without limit as one approaches the critical solution''. For branes that reach the 
center, 66' = at the center and hence TZdt scales as 56~'^. For solutions ending away from 
the center, 56' diverges, and the curvature at the endpoint scales as p~^. Under a scaling 
transformation, both p and 56 scale as p, so the curvature in either case will scale as jU~^. 

When the scalar curvature of the induced metric becomes of order of the string scale, 
a'~^ = stringy corrections to the DBI action will no longer be negligible, and hence our 
analysis of the scaling behavior (based entirely on the DBI action) will cease to be valid. 
In other words, higher order corrections will become important when p~'^ ~ These 
corrections are due to excited open string modes whose mass naively goes as £~^, but whose 
mass is reduced by a power of p when the endpoints of the string are in the high-curvature 
region. To illustrate this issue graphically, we plot ln[— T^/jy] in figure 12, evaluated at the 
point of closest approach to the center, as a function of the asymptotic coefficient 0(o)- 

Using a'~^ = ^ ^ ^ ^ ^q^^^ _ 6*^^^^)^/^ = {m- m*)^/^, the condition p~'^ ~ 
becomes (m — m*) ~ A^^^^. The properties of the boundary theory are thus governed by 
the critical exponents we calculated above, but only in the range of masses \~^^^/R?j <C 



'^For the D7-brane in x slicing (with its first order transition), the scalar curvature also diverges for 
the critical solution, but in this case the critical solution lies on an unphysical (infinitely unstable) branch. 



-27- 



(m — m*) <ti l/i?3. (Once again, the same relations are true for m < m* , with m — m* 
replaced by m* — m, etc.) In particular, in this window the three-point function of the zero- 
momentum Om operator grows as (m — m*)~'^^^. This power law growth will eventually be 
cut off at the lower end of the scaling window. As A is a free parameter, however, we are free 
to consider the regime where this scaling window extends over arbitrarily many decades. 



5. Conclusion 



The scaling symmetry of near-center or near-horizon probe Dp-brane solutions in global AdS 
has allowed us to determine when phase transitions associated with fundamental-representation 
fields coupled to A^ = 4 SYM theory will be first order or continuous. When they are continu- 
ous, we find that the approach to criticality is governed by non-trivial critical exponents which 
can be calculated analytically using the supergravity description. We emphasize that the 
phase transitions we have found are finite- volume, large- A'^c, and large-A effects. The continu- 
ous transition we find in the D7-brane case can be interpreted as a meson binding/unbinding 
transition similar to the finite-temperature case analyzed in previous studies. 

A finite volume in the flavored SYM theory can give rise to interesting new effects in 
the phase diagram, such as the appearance of our continuous phase transition. We have only 
explored a small part of the multi-dimensional phase diagram of this theory. Probe Dp-brane 
techniques can be used to study systems at finite density [24, 36-39] or in background electric 
and magnetic fields [19,20,40-44]. We expect a rich phase structure to emerge for these 
systems when confined to finite volume. 
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